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ABSTRACT 

In this paper, a general model of wireless channels is established based on the physics of wave propagation. Then 
the problems of inverse scattering and channel prediction are formulated as nonlinear filtering problems. The 
solutions to the nonlinear filtering problems are given in the form of dynamic evolution equations of the estimated 
quantities. Finally, examples are provided to illustrate the practical applications of the proposed theory. 
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1. INTRODUCTION 

In modern engineering, the problem of modeling wave propagation in a random field of scattering objects and 
of the associated signal processing, has become more and more important, with wide applications ranging from 
wireless communications, to radar or sonar object detection, to medical imaging using microwave or infrared 
light. In all of the applications, the first thing is obviously to understand and model the wave propagation 
in random fields. In the object detection and imaging applications, the signal processing is usually an inverse 
scattering problem, that is to find out the configuration and dynamics of the scattering objects in the random 
field, based on the data of received signals and the known transmitted signal. When applied to the field of 
wireless communications, for instance, channel estimation and prediction for dynamic power control or adaptive 
modulation, the problem goes one step deeper and becomes filtering and predicting the dynamics of the channel 
response, given some observed data of the past channel response. A powerful model of wave propagation will be 
established based upon the physics of wave scattering, by which the channel response may be computed explicitly 
as a functional of the random field of scattering objects. Then, it becomes natural to approach the problem of 
estimating and predicting the channel response in two steps: solving the inverse scattering problem first, then 
calculating the channel response using the solved dynamics of the scattering objects. 

Such a natural two-step approach seems to have been overlooked in the existing literature of communications. 
In fact, the traditional approaches try to avoid the inverse scattering problem. They often start directly from 
some presumed stochastic models of the channel response, without looking into the actual random field of 
scattering objects at all. In the Rayleigh fading model that has been widely used since the early days of wireless 
communications,^ the physical scattering environment is absent from the consideration, while the received signal 
is assumed to be a sum of many differently delayed and weighted versions of the transmitted signal, where the 
weighting coefficients are complex Gaussian and independent, the delay times modulo l//c are also independent 
and uniformly distributed in [0, l//c), with fc being the center frequency of the carrier.^ When there is one 
strong, dominating component, e.g. from line-of-sight propagation, the Rayleigh model is slightly modified into 
the Ricean fading model. Both models are based on the assumption of a large number of scattering objects, 
which may not always be the case. In case an adaptive mechanism is employed to control the transmitter power 
or vary the transmission bit rate, the channel estimation and prediction circuit has to be very fast to follow the 
dynamics of the channel response.^ The previous channel models have become inadequate in millimeter wave 
mobile communications, where the channel is fast-fading, namely, the bit duration is short comparing to the 
dynamics of channel response.^ In fact, the inadequacy of the old models may be held responsible for a general 
opinion among many wireless engineers that in many fast-fading systems, it might be impossible to estimate the 
channel state information (CSI) so to adapt the system accordingly.^ Such a pessimistic view has even inspired a 
lot of work concerning the capacity and coding schemes of future high speed and fast-fading wireless channels.'^' ^ 
A generally used estimation of the fading speed is based upon the assumption that the signal fading experienced 
by two receivers become uncorrelated when the separation between the receivers is significantly larger than the 
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carrier wavelength.^ However, in the Rayleigh or Ricean model, the signal fading is governed by a Gaussian law 
of distribution, which implies the equivalence between uncorrelation and independence. So even in principle, the 
signal fading seen in one region does not contain any information for the signal fading in another region just a 
few wavelengths away. This imposes a serious limitation upon the estimation of CSI. For instance, consider a 
wireless link with the carrier frequency around 1.9 GHz, so the wavelength is about 0.158 meter. It takes merely 
8 milliseconds for a mobile with speed 20 m/s to travel through one wavelength. It seems, according to the 
Rayleigh and Ricean models, that an estimator for the CSI has to accomplish its task well within a millisecond, 
and the obtained CSI can be regarded as reliable for no longer than 8 milliseconds. 

We think that the above limit on the CSI estimation of a fast-fading channel is over pessimistic. Our point 
of view is that for any wireless link in a given physical environment, the channel response is the result of the 
superposition of waves propagating via multiple paths, where each signal path is a geometric route that starts 
from the transmitter, possibly hits one or more scattering objects, and finally reaches the transmitter. The 
physical environment that contains the scattering objects is called a scattering environment. The response of a 
wireless channel is actually deterministic, conditioned upon the scattering environment and the locations of the 
transmitter and receiver inside the scattering environment. The randomness of wireless channels arises only from 
the variety of scattering environments. In a wireless channel model based on the physics of wave propagation in a 
scattering environment, a natural choice for the "channel state" would be the state of the scattering environment, 
namely the positions and velocities of the scattering objects as well as their physical effects on the wave. So the 
estimation of CSI boils down to estimating the positions, velocities, and attenuation coefficients of the scattering 
objects. Such channel state should remain largely time-invariant over a much longer time period. As an example, 
consider a scattering environment in which the distance among the scattering objects, the transmitter, and the 
receiver is at least several tens of meters, so that the mobile receiver or transmitter experiences the same scattering 
environment even after it travels over 10 meters. With the 20 m/s mobile speed, it takes 0.5 second to cover 
that distance, which is more than 60 wavelengths for a carrier frequency around 1.9 GHz. 

There have been previous works adopting a deterministic model of signal fading, ^"^^ but in most of them 
the effect of the scattering environment is represented by a siipcrposition of plane waves in proximity to the 
receiver, whose interference pattern is sampled by the receiver. This representation only works for a frequency- 
nonselective channel,^ where the frequency bandwidth of the signal is too narrow to resolve the delay times 
of different paths. When the signal bandwidth is sufficiently large to resolve the path delays, the plane wave 
representation breaks down, a more advanced model that captures both the direction of arrivals and the path 
delays is necessary. Moreover, the previous works^"^° estimate the CSI using various spectral estimation tech- 
niques, which do not guarantee the optimality of the solutions. In this paper, we shall first establish a complete 
model of wireless channels based on the physics of wave propagation. Then the problems of inverse scattering and 
channel prediction are formulated as nonlinear filtering problems. The solution to the nonlinear filtering problem 
will be given in the form of a dynamic evolution equation (the filtering equation) for the estimated quantities, 
just like the Kalman filtering equation for the linear filtering problem. In most practical applications, the 
filter equation may be implemented, or at least be approximated, by a state-possessing machine driven by an 
innovation input. In our particular applications, it happens that both the state of the scattering environment 
and the observation process are governed by linear differential equations. Note that even though the system 
equations are linear, the Kalman filter may still be sub-optimal, unless the initial probability law of the system 
state is Gaussian. It is proven mathematically that the nonlinear filter produces the optimal estimation for the 
interested quantity. Even the nonlinear filter may not be easily implemented in practice, its solution provides a 
benchmark for evaluating the performance of other sub-optimal estimators. The theory of nonlinear filtering has 
been well established, many textbooks exist, covering its fundamental principles and applications, especially to 
the field of control engineering, ^^"^^ however it has not made much appcaranc;e in the field of c;omnimiication 
engineering. Apart from an attempt to solve the channel prediction problem for fast-fading wireless channels, 
the present paper is intended to introduce the nonlinear filtering theory as a powerful tool for signal processing 
and system optimization in the field of communication engineering. 

2. RANDOM SCATTERING FIELDS AND WAVE PROPAGATION THEREIN 

Throughout the present paper, a sufficiently large probability space (O, P) is assumed in which all random 
variables under consideration are defined. In the wave-involved applications mentioned in the beginning of the 
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introduction, there is often a random scattering field (RSF) consisting of scattering objects, due to which the 
signal wave from a transmitter to a receiver undergoes multiple scattering and propagates along various geometric 
paths associated with different power losses, different time delays and direction of arrivals, as pictured in Fig.l. 
The spatial domain £ C H'? accommodating the scattering objects together with the signal transmitter and the 
receiver is called the scattering environment, where g G IN is the dimension of the space under consideration, 
usually g = 2 or 3. In general, an RSF, denoted by iS'(i;r,v), may contain spatially extensive scattering objects, 
which could be continuously distributed in the real space and the velocity space, namely the configuration space. 
Fig.l only shows a special RSF with discrete scattering objects, which may be the most adopted model. A 
discrete RSF is often modeled by a compound point process^^' "'^^ in the configuration space £ x IR*^, 

N[t) 

^(t;r,v) = ^5fc,5(r-rfe)(5(v-v,), (1) 

fc=3 

where r — (ri, • • • , r^)' G £ and v = {vi, - ■ ■ ,Vq)' G M'^ denote* the position and the velocity respectively, and 
hence {(rfe, Vfe)}^^'' are configuration coordinates for the scattering objects, which are the realization points of 
a space-time point process in [0,cxd) x £ x IR'^. Here N{t) is the total number of points occurred up to time t 
anywhere in £ x IR'^. The scattering objects are indexed by the integers A: > 3, the indices k = 1,2 are reserved for 
the transmitter and the receiver respectively. For each fc > 3, S**; is a complex- valued random variable associated 
with the kth realization point, which represents the scattering response of the corresponding scattering object 
to the signal. Although the discrete RSF model is chosen here to simplify the discussion, it should be pointed 
out that, a more general RSF model might, and should be used when necessary, for example, when describing a 
scattering environment with continuously distributed scattering objects. In fact, all mathematical formulae with 
iS'(i;r,v) involved in the present paper arc formulated in such general forms that they hold for any !)'{£ x IR')- 
valued process S{t;r,v), where 'D'{£ x IR'^) denotes the Hilbert space of the complex- valued, linear functionals 
of functions in £ x H'', which are usually called distributions or generalized functions. 
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Figure 1. A typical random scattering field, in which the signal wave propagates from a transmitter to a receiver along 
various geometric paths as a result of the interaction with the scattering objects. 

The RSF imprints itself on the channel response between a wave transmitter and a signal receiver, all of 
which lie in the scattering environment £. In general, both the transmitter and the receiver could be extensive 
in space, as being the case when antenna arrays are employed. Therefore we use a;(t;ri,vi) and y(t;r2,V2) to 
denote the spatially extensive signal amplitudes of the transmitter and the receiver respectively, both of which 
are complex- valued and square integrable in the Lebesgue sense, namely, x,y G L^([0, cxd) x £ x IR'',C), where 
C is the set of complex numbers. Recalling the theory of wave scattering which says that the otherwise freely 

*In this paper, A' denotes the transpose of A, for any vector or matrix A. 
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propagating signal wave from a source is perturbed by the scattering objects, and the actual wave in space is the 
superposition of the free-running wave and all of the scattered waves. The Huygens-Presnel principle^^' further 
states that the scattered wave from each scattering object, often called the second wave, is again a free-running 
wave coming out of the object, with an amplitude proportional to that of the incident wave, and the proportional 
coefficient is called the scattering response of the object. In general, it is rather complicated to calculate the 
actual wave in an exact manner, which involves the solution of an integral equation or a series expansion using 
the Feynman diagram. The complication has to do with the multiple scattering effects. However, when 
the scattering objects are sparsely distributed in space, their effect on the signal wave may be easily calculated 
using the perturbation theory. In particular, the first-order perturbation theory approximates the actual wave by 
superposing the free-running wave and the second waves excited by it at the scattering objects. For an obvious 
reason, the first-order perturbation is also called the single-bounce approximation. When considering only the 
contribution of the second waves, the single-bounce approximation results in an input-output relation of the 
spatially extensive signals a;(t;ri,vi) and y(t;r2,V2) in their low-pass equivalent form,^ 

Nit) 

2/(t;r2,V2) = ^ / a;(t - Tife - rfe2;ri,Vi)/(cTife)S'fc/(cTfe2)exp(i27r7ife2t)(iridvi, (2) 

wherc rife and Tk2 arc the time delays associated with the signal propagation from the transmitter to the fcth 
object and from the kth object to the receiver respectively, which satisfy the following equations, 

|ri -rfc(^-Tfc2)| |rfcft-Tfe2) -r2| w , . „ 

rife = , rfe2 = , V K > 6. [6) 

c c 

Also in the above equations, c is the speed of the wave, 1{ct) is the propagation loss that depends only on the 
propagation distance cr, and 

def /o[vi - Vfe(i-Tfe2)] • [ri -rfe(t-Tfe2)] , fo[Vk{t - Tk2) - V2] ■ [rk{t - Tk2) - r2] ... 
7U2 = j 77 71 \ i 77 7 \ ) \V 

c|ri - rfe(f - Tfe2)| c|rfe(i - rfe2) - r2| 

is the Doppler frequency shift experienced by the signal with the central carrier frequency /o. Note that the time 
delays {Tfe2}fe>3 are implicitly defined in equation (3), although an explicit formula may be derived when the 
kinetic motions of the scattering objects are given explicitly. For example, if the velocity difference Vfe(t) — Vfe(t — 
Tfe2) is neglected for all fc > 3, that is to ignore the velocity variation during the time of the wave propagating 
from the scattering object to the receiver, then the following equation holds, 

crfe2 = \rk{t) - Vfc(t)rfe2 - r2|, V fc > 3, (5) 

which is easily solved to give an explicit formula, 



Then the time delays {Tik}k>3 arc explicitly computed using equations (3) and (6). Similarly, equation (4) could 
be simplified if all {vfe}fe>3 are treated as constants during the interval [t — Tk2,t], 

def /o[vi - Vk{t)] ■ [ri - rfe(t) + Vfe(t)rfe2] , /o[vfe(t) - V2] ■ [rk{t) - Vfc(<)rfe2 - r2] w , . „ /'7\ 

= \ 77^, 7T\ 1 ^ 1 — 777 777 , > V > 3. (7) 

c|ri - rfe(t) + Vfc(t)Tfc2| c\rk(t) - Vk(t)Tk2 - r2| 

Neglecting the object velocity variation during the time of wave propagation is often an excellent approximation 
to make, especially when dealing with electromagnetic waves, which travel at the ultimate speed in nature, 
c = 3.0 X 10** m/s. Even for ten kilometers of propagation distance, the time delay is merely 1/30 millisecond, 
during which no macroscopic scattering object can change its speed for much. Furthermore, the highest speed 
of most macroscopic objects is no more than a millionth of the light speed. Consequently, equation (6) may be 
simplified, 

^^^^|rfe(0-r2|_M0-r2]-Vfeft)^^^^3^ (8) 
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by neglecting terms of the second or higher powers of Vk(t)/c, V fc > 3. Or even simpler, while still an extremely 
good approximation for the propagation of electromagnetic wave, the time delays may be calculated as, 

Tlfe = , Tk2 = , V K>6. (9) 

C C 

To write the input-output relation in a more compact form, we define 



'r2(r,v;r2) = ^2 _ |v|2 ' V A; > 3, (10) 

, X def ki -r + vr2(r,v;r2)| „ 
Ti(ri;r,v;r2) = ^, V A; > 3, (11) 

^tr- .r.r...r. ^ def /o(vi " v) ' [ri - r + VT2 (r, v; r2)] , /o(v - V2) ■ [r vTa (r, v; r2) - r2] 

7iri,Vi;r,v;r2,V2) = ■ ■ ^ 1 , r ■ , (12) 

c|ri - r + vT2(r,v;r2)| c|r - vr2(r, v; r2) -r2| 

L(ri, vi; r, v: r2, V2) =^ /[en (ri ; r, v; r2)]/[cT2(r, v; r2)] exp[i27r7(ri, Vi; r, v; r2, V2)], (13) 
then we use equation (1), together with equations (6), and (7) to rewrite equation (2) as, 

2/(t;r2,V2) = / a;[t-Ti(ri;r,v;r2) -r2(r,v;r2);ri,vi] 

X L(ri,vi;r,v;r2,V2)S'[t - T2(r,v;r2);r,v]dridvidrdv. (14) 

This equation shall serve as the fundamental input-output relation describing the propagation of spatially exten- 
sive signals in an RSF, up to the single-bounce approximation. It is evident that the received signal y{t; r2, V2) 
is a bilinear fimctional of a;(t;ri,vi) and S{t;r,v), for any fixed (t;r2,V2). More specifically, j/(t;r2,V2) is a 
linear functional of x{t; ri, vi), if L(ri, vi; r, v; r2, V2)5(t; r, v) is regarded as an integration kernel; or y{t; r2, V2) 
could be viewed as a linear functional of S{t; r, v), with x{t; ri, Vi)L(ri, vi; r, v; r2, V2) serving as an integration 
kernel. 

3. THE DYNAMICS OF RSF 

Let us consider the case in which the set of scattering objects is fixed in time and the scattering response of each 
object is time- independent. Nevertheless, the RSF is still time-varying when the objects are subject to kinetic 
motion. We shall employ an RSF model so general that the scattering objects are subject to accelerations 
that depend upon the position. Specifically, the motion of the fcth scattering object. A; > 3, is described by a 
differential equation, 

dTk = Vkdt, (15) 
dvk = f{rk)dt, (16) 

where /(r) is an M'- valued function representing the deterministic acceleration exerted upon the scattering 

objects. Since the total number of scattering objects within the environment £ is assumed to be fixed in time, 
it is a random variable A'' = N{lo). As mentioned before, all random variables and processes are defined on the 
"grand" probability space {fl,J^,P). Let 

N 

S{t; r, v) ^^^("^ - ^'^)^(^ - ^fe)' (1^) 

fe=3 

then 5(t;r, v) is a space-time compound point process^^'^* in the space [0, 00) x £ x R^. Let U denotes the 

Hilbert space of all generalized functions in £" x IR«, namely U = V\£ x W), then S{t] r, v) may also be viewed 
as a ?7- valued random process indexed by t e [0,oo). Next we attempt to derive a partial differential equation 
for 5(4; r,v) using equations (15), (16) and (17). 
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While explicitly parameterized by r and v, the RSF S{t; r, v) of equation (17) is implicitly a function of the 
processes {(r/s, Vfe)}^3, described by the differential equations (15) and (16). By differentiating equation (17), 
it is obtained that, 



dS{t;r,v) = Sk5{\- - Vk)d5{r - rfc) + Sk6{r - rk)d6{-v - Vfc) 

k k 

^^^(v - Vfc) — -Vkdt + SkS{r - Ffe) — ^— 



^SkS{v -Vk) 



drk ^ dvk 

k 

-Vkdt - ^Skd[r - Tk)- 
k 



f{vk)dt 



dr 



dv 



f{rk)dt 



■ E SkSiv - v.)^fc^vdt - E SkSir - rk)^^^^^mdt, 



for all A; > 3, where in the last step, two identities have been used, 

<5(v-Vfc)vfe = (5(v-Vfe)v, VA;>3, 
5{v-Vk)f{vk) = (5(r-rfe)/(r), Vfc>3. 

The partial differential operator d/dr is defined as a row vector. 



dv 



(18) 



(19) 
(20) 



def 

dv 



d_ d_ 

5 ' * * ? 



dr„ 



(21) 



and the operator d/dv is similarly defined. Note that /(r) is valued as a g-dimensional column vector. When a 
row vector is followed by a column vector, it is understood that the scalar product in the space IR' takes place. 
Noticing two more identities. 



dS{t;v,v) ST o xi .dS{r-rk 

= 2^ ^^^^(V - Vfc) 



r '^-^ or 

k 

dS{t;r,v) \- q r/^ ^ ^ dS{v-Vk) 



(22) 
(23) 



one may easily simplify equation (18) into a partial differential equation governing the RSF, 

,r,,. ^ dS{t;r,v) dS{t;r,v) 

dS{t;r,v) = '-vdt ^— -f{r)dt. 

or ov 

Define a linear operator A: U ^ U, such that for any ^ 



^ def OS, 9^ 
Ae(r,v) = --v--/(r), 



(24) 



(25) 



then equation (24), together with a proper initial condition, is immediately recognized as an infinite dimensional 
linear system^^' in the Hilbert space U, 



dS{t) = AS{t)dt, (26) 
S{0) = So, (27) 

where is an (.F, ;B(C/))-measurable random variable, B{U) is a Borel cr-algebra of subsets of U. It is to be 
stressed again that the dynamic equations (26) and (27) are applicable to any RSF S{t) G C/ = V'lS x IR'), 
although the derivation of the equations is exemplified by the special case of discrete RSF for simplicity. 
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4. SIGNAL DETECTION 

As discussed in section 2, when a spatially extensive signal x{t; ri , vi) is transmitted, a spatially extensive receiver 
will catch a signal j/(t;r2,V2), which may be degraded by a white Gaussian noise. Let 



i) = y{ 

Jo 



r(t;r2,V2) = / y{s;r2,V2)ds + Z{t;r2,V2), (28) 



where Z{t; r2, V2) is an iJ-valued Wiener process representing the observation noise, with H =^ L'^{£ x M'), that 
is, the space of square Lebesgue integrable functions in f x IR"^. It is also assumed that the covariancc operator 
of Z{t), denoted by R, is uniformly positive definite. Using equation (14), with the simplified calculation of delay 
times in (9), we can write. 



/ |ri -r| + |r-r2| , 
x[t ;ri,vi ) L(ri,vi;r,v;r2,V2) 

(£-xR9)2 V 

sit — J -^;r, V ) dridvidrdv 



dY{t;r2,V2) = 

dt + dZ{t;r2,V2), (29) 

F(0) = 0, (30) 

which is the signal detection equation, or called the observation equation. With the observation data {Y{s)}o<s<t, 
two estimation problems emerge, depending upon the actual application goal: 

1) estimating and predicting S{t), when {a;(s)}o<s<t is known; 

2) estimating x{t), when {S{t)}o<s<t is known. 

The first problem is related to various applications of inverse scattering, such as radar object detection, channel 
estimation and prediction for wireless communications; while the second problem is specific to signal detection, 
especially to the optimal design of a receiver in communications. In the present paper, we arc interested in the 
first estimation problem, in which the signal x(t; ri, vi) is regarded as a known quantity. Define a time-dependent 
linear operator G{t) : U ^ H, such that for any time- varying function ^{t) G f/, V t > 0, 

[ x(t-— ^hL^ — ^;ri,vi ] L(ri,vi;r,v;r2,V2)^ U - -^^^ — ^; r, v ] dridvidrdv, (31) 

■/(£xRo)2 V C J \ C J 

then the observation equations (29) and (30) are of the form, 

dY{t) = G{t)S{t)dt + dZ{t), (32) 
Y{0) = 0, (33) 

where Y{t) and Z{t) are H-valued random processes, with H = L^{£ x IR'), while G{t) is a time-dependent 
linear operator in Chs{U, H), with the time-dependence originated from the known and time- varying signal x{t), 
^Hs{U,H) being the space of all Hilbert-Schmidt operators from U to H. 

The system equations (26), (27) and the observation equations (32), (33) constitute a nonlinear filtering 
problem, that is to compute the conditional probability, 

n(i)(V(i)) =E[^{t,SmY{s),0< s<t], (34) 

where 1^(1,^) G C^'^([0, 00) x t/) is a given test function, which is boimdcd and has continuous first and second 
order Frechet derivatives with respect to t and ^ respectively. The essence of the nonlinear filtering problem is 
to obtain the optimal estimation of the interested quantity 'ip{t, S{t)), based on the observation data up to time 
t, {Y{s)}o<s<t- In most practical applications, it is suflacient to obtain a dynamic evolution equation for the 
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estimated ip{t,S(t)), usually in the form a differential equation driven by the observation data, which can be 
implemented, or simulated, by a state-possessing machine or circuit. Accordingly, a dynamic evolution equation 
for the estimated quantity may be regarded as a solution to the nonlinear filtering problem. Thanks to the greatly 
advanced theories of nonlinear filtering and of finite and infinite dimensional systems, such dynamic evolution 
equations have been well established as solutions to the nonlinear filtering problems. 

5. CHANNEL ESTIMATION AND PREDICTION 

For the general theory of nonlinear filtering, the readers are referred to the excellent textbooks^^"^^ and references 
therein, especially the original papers. Here we shall formulate the problem in a less general form, but to suit 
our applications better. We are interested in an inverse scattering problem, which is recast as, given, 

dS{t) = AS{t)dt, S{0) = So, (35) 

dY{t) = G{t)S{t)dt + dZ{t), F(0) = 0, (36) 

to compute, 

mm)) = E[^(i, smY{s), o<s<t], (37) 

where Z(t) is an ff-valued Wiener process, with covariancc operator i?, which is uniformly positive definite. Sq 
is an (JF, S({7))-measurable random variable with probability measure Hg. Z{t) and Sq are independent. The 
operators A and G{t) are defined in (25) and (31) respectively. Although we have successfully formulated the 
general inverse scattering problem as a nonlinear filtering problem for an infinite dimensional linear system, we 
shall not present the general solution to the problem, in order to limit the scope and length of the discussion. 
The general solution involves much elaborate notations and concepts in functional analysis, as well as somewhat 
tedious arguments of convergence, which does not seem to fit the purpose of the present paper. Interested 
readers, however, are referred to the literature.^'"' 

What we shall focus here, is a special case of discrete RSF with the accelerations of the scattering objects 
being neglected, namely, all objects are moving with constant velocities. The goal of inverse scattering is to 
find out the initial position and velocity v^, as well as the attenuation coefficient Sk to the wave, for each 
scattering object k, 3 < k < N, based upon the received signal {Y{s)}o<s<t and the known transmitted 
signal {x{s)}o<s<t- However, the total number of scattering objects TV is usually random, and to pin down 
the configuration coordinates {(ffc, Vfc)}j^2 from scratch is highly difficult. A more practical approach is to 
discretize the initial configuration space into small regions and label each of them by an integer n, 1 < n < M, 
then to assign a complex-valued random variable S„ as the scattering response to the center of the nth region, 
which is coordinated by (r„,v„). The idea is, when the initial configuration space is divided into suSiciently 
small regions, the scattering objects within any region n, 1 < n < M, become non-resolvable, their scattering 
responses superposed together is represented by a single scattering point sitting at the center of the region with 
a scattering response S„. The result of the discretization is to remove the coordinates of the scattering objects 
in the initial configuration space from the set of the unknowns, but to let the complex random variables {Sn}n=i 
bear all the information about the scattering environment. The probability law of {Sn}n=i may be obtained 
from some experimental data, or calculated from an empirical model of the distribution of scattering objects in 
an actual geographical environment.^^' For instance, if the distribution of objects is quite sparce in an area, 
then the probability of S„ = is high, for any 1 < n < M. Let 

S = (Si, S2, • • • , Sm)', (38) 

then S is a random variable in C^^ , whose probability law is known and denoted by Hq just as before. Since the 
scattering response of each object is assumed time-invariant, the RSF S{t), now C^-valued, is time-independent, 

dS{t) = 0, S{0) = S. (39) 

The transmitter and receiver are also discrete in the configuration space, which consist / and J antenna elements 
located at {{rii,yii)}j^i and {(%, V2j)}^^i respectively. As a consequence, the transmitted signal x{t) is a 
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C'^-valued signal, while the Wiener noise Z{t) at the receiver and the received signal Y{t) are C'^-valued random 
processes, 

x{t) {xi{t),X2{t),---,xj{t)y, (40) 

Yit) {Y,{t),Y2{t),---,Yj{t)y, (41) 

Z{t) {Z,{t),Z,(t),---,Zj{t)y. (42) 

Again let R{t) denote the covariance matrix of the Wiener process Z{t), which is uniformly positive definite. 
Define a J x M matrix G{t) = {gjnit)}i<j<js<n<M, such that, 

9jn{t) = E ^ (* - ''"'"""'^ rii, vii) i(rH, vu; r„, v„; rsj, V2,), V 1 < j < J, V 1 < n < M, (43) 

where the function L is defined in (13), then the discrete version of the observation equation (29) is formulated 
as, 

dY{t) = G{t)Sdt + dZ{t), y(0) = 0. (44) 

With the trivial system equation (39) and the simple observation equation (44), the inverse scattering problem 
is to compute, 

n(t)(^(t)) =E[7/.(t,5(t))|r(.s),0 < s < t]='E[yj{t,S)\Yis),0 < s < t], (45) 

where il^{t,S,) G C^'^([0, cx)) x C^,C) is a given test function, which is bounded and has continuous first and 
second order derivatives with respect to t and ^ respectively. For the convenience in presenting the following 
filtering equations, let B : [0, oo) x C'^ C'^ he a function such that, 

B{t, G(i)C, V i > 0, V ^ e C^, (46) 

B'{t, ^'G'{t), V t > 0, V C e C^. (47) 

It is a trivial matter to verify the conditions in order to apply the Kushner-Stratonovitch equation,^^' 

dnmrn) = u{t)(^)dt+[ii{t){B'it)m)-mim)m{B'{m 



x R-'^{t)[dY{t) -n{t){B{t))dt\, almost smely, \/ t> 0, (48) 

n(o)(v(o)) = no(V(o))=/ ^(o,Orfno(0. (49) 

Another well known solution to the nonlinear filtering problem is the Zakai equation, ^'^ which governs the 
so-called unnormalized conditional probability p(i)(^(t)), 

dp{t){ip{t)) = p{t) ( ^ ) dt + p{t){B'{t)7p{t))R-\t)dY{t), almost surely, V t > 0, (50) 



p{o){m) = Mm)= f modMo- (51) 

The relation between p{t) and Il{t) is given by, 

n(i)(V') = ^ test function ip, (52) 

where p{t){l) satisfies, 

dp{t){l) = p{t){B'{t))R-\t)dY{t), almost smely, y t>0, (53) 

p(0)(l) = 1. (54) 
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As a simpler example, consider a flat-fading channel from a single transmitter element to a single receiver 
element, where the frequency bandwidth of the channel is sufficiently narrow that the time delays among the 
signal paths are insignificant, however, the channel response experiences a fast oscillation due to the motion of the 
receiver, or the transmitter, or the scattering objects. The goal is to estimate the channel state and to predict 
the channel response, based on the known transmitted signal and the data of received signal. This problem 
has been tackled before with the help of various spectral estimation techniques, ^^^"^ which do not guarantee 
the optimality of the solutions. It is therefore of great interest to formulate the problem using the language of 
nonlinear filtering, and to see what kind of optimal solution may be obtained. The observation equation, namely, 
the channel input-output relation, has a very simple low-pass equivalent form,^ 

dY{t) = x{t)h{t)dt + dZ{t), y(0) = 0, (55) 

where x{t) and Y{t) arc the transmitted signal and the integrated version of the received signal respectively, 
and h{t) is the time-varying channel response, all C-valued, while Z{t) is a C-valued Wiener process with a 
uniformly positive definite covariance function R{t). The channel response h{t) is a superposition of many 
Doppler components, 

K 

Kt) = ^^'^k exp{i2wfkt), - /max < fl < h < ■ ■ ■ < Ik < /max, (56) 

where /max is the maximum Doppler frequency shift, and {fk\^=\ constitutes a partition of the frequency band 

[—/max, /max], and Ofc G C is thc constant amplitude of the fcth Doppler component, \ < k < K . It is assumed 
that the partition {Jk}k=i sufficiently fine to justify the discrete approximation to the Doppler shifts. Define 
a function F : [0, oo) x — > C, such that, 

K 

F{t,0 = x(t) ^exp(i27r/,i)a, V C = (a, 6, • • • ,Cif)' e C^. (57) 

k=l 

The nonlinear filtering problem for channel estimation and prediction is formulated as, given, 

a '= (oi, 02, • • • , qk)' G , with probability law XIq, (58) 

dY{t) = F{t,a)dt + dZ{t), ^(0) = 0, (59) 

to compute, 

n(i)(V(t)) = E[V(t, a)\Y{s), 0<s<t], (60) 

where e C^'^([0,oo) x C^,C) is a test function. Again, the solution is given by the Kushner-Stratonovitch 
equation, 

du{tmt)) = u{t)(^^yt+[u{t)iF{t)m)-mim)m{Fm 

X R-\t)[dY (t) -n{t){F{t))dt], gilmostsmely, V t>0, (61) 

n(0)(^(0)) = no(V'(0))= / ^iO,OdILoiO, (62) 
or by the Zakai equation, 

dp{t){il>{t)) = p{t) (^^^ dt + p{t){F{t)^{t))R-'^{t)dY{t), almostsmely, y t> 0, (63) 

pmm) = Mm)= [ modMO, m 

then the probability Il{t){il){t)) is calculated as, 
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6. CONCLUSION 



It has been argued that the input-output response of a wireless channel in a given scattering environment is 
actually deterministic in nature, conditioned on the positions and velocities of the scattering objects, as well 
as the locations and velocities of the transmitter and receiver in the environment. Using a physical model of 
wave propagation, the channel response is calculated explicitly as a functional of the random scattering field. 
Consequently, the best choice of the channel state information for wireless channels may be the state of the 
random scattering field, which can remain unchanged for a sufficiently long time, even for the "fast-fading 
channels" in the conventional sense. As far as the problem of channel estimation and prediction is concerned, 
there seems to be no fundamental limit imposed by the physics of wave scattering, no matter how short is the 
carrier wavelength and how fast the mobiles move as long as the inter-object separation is much larger than the 
distance that any mobile can travel within the time duration of interest. Only practical limits may rise from the 
speed and complexity of the actual signal processing circuits. The problems of inverse scattering, formulated as 
nonlinear filtering problems, are solved in terms of dynamical evolution equations for the estimated quantities. 
However, the filtering equations cannot be implemented with little difficulty using today's signal processing 
circuits. Both better numerical algorithms and more advanced electronics are desired, in order to take fully the 
advantage of the nonlinear filters in practical applications. 
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